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ABSTRACT 

Using the energy variational method of magneto-sohd-mechanical theory of 
a perfectly conducting elastic medium threaded by magnetic field, the frequency 
spectrum of Lorentz-force-driven global torsional nodeless vibrations of a neutron 
star with Ferraro's form of axisymmetric poloidal nonhomogeneous internal and 
dipole-like external magnetic field is obtained and compared with that for this 
toroidal Alfven mode in a neutron star with homogeneous internal and dipolar 
external magnetic field. The relevance of considered asteroseismic models to 
quasi-periodic oscillations of the X-ray flux during the ultra powerful outbursts 
of SGR 1806-20 and SGR 1900+14 is discussed. 
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1. Introduction 

In the context of recent discovery of quasi-periodic oscillations (QPOs) in the X-ray 
luminosity during the giant flare of SGR 1806-20 and SGR 1900+14 that were interpreted 
as being produced by torsional vibrations of quaking magnetars (Israel et al. 2005, Watts 
& Strohmayer 2006), in (Bastrukov et al. 2009a, 2009b) several scenarios of the post-quake 
vibrational relaxation of a neutron star model with uniform internal and dipolar external 
magnetic field 

Br = Bcos9, Bg = -Bsm9, B^ = 0, r<R (1) 
RV . ^ B fR^ ^ 



Br = B { - j cosO, Be = -—ij] sine, B^ = 0, r> R (2) 
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have been studied on the basis of equations of Newtonian magneto-sohd-mechanics 

pu=-[5jxB], 6} = —[Vx6B], 5B(r, t) = V x [u x B]. (3) 

c An 

These equations describe the Lorentz-force-driven non-compressional {6p = — pV ■ u = 0) 
fluctuations of star matter about axis of above fossil magnetic field B in terms of fluctuating 
material displacements u (the basic dynamical variable of solid mechanics) and the magnetic 
field (5B. It is implied that elastic stellar material is of an extremely high electrical conduc- 
tivitjEI. The chief argument for interpreting QPOs during the outbursts of above mentioned 
SGRs (detected on descending branch of the giant flare light-curves) as caused by torsional 
vibrations of a solid star driven by restoring force of magnetic field stresses is that it is an 
ultra strong magnetic field frozen-in the entire volume of magnetars serves as the main en- 
ergy source and promoter of their X-ray bursting seismic activity (e.g., Woods & Thompson 
2006, Mereghetti 2008) and, also, bearing in mind the fact that the very notion of torsional 
vibrations has come into theoretical seismology from the solid-mechanical theory of shear 
vibrations of an elastic sphere (e.g. Lapwood, & Usami 1981, Lay & Wallace 1995, Aki & 
Richards 2002, Stein & Wyssesson 2003). It worthy noting that theoretical investigations of 
non-radial torsional Alfven oscillations of a fluid star in its own homogeneous magnetic field 
have a long story that was started, to the best of our knowledge, in works of Jensen (1955) 
and Plumpton & Ferraro (1955). Remarkably that in the latter paper, by emphasizing the 
basic discovery of Alfven that perfectly conducting fluid threaded by magnetic field behaves 
like anisotropic elastic medium capable of transmitting mechanical disturbance by transverse 
hydromagnetic waves (e.g. Falthammar 2007), it is argued that eigenfrequency problem of 
such vibrations can be tackled on the basis of equation 

pii = ^[V X [V X [u X B]]] X B, V-u = 0, V-B = (4) 

An 

which, as is evident, follows form equations ([3]). It is clear from this last equation that the 
frequencies of Alfven oscillations must substantially depend on both geometrical configura- 
tion of internal equilibrium magnetic field B and analytic form of oscillating field of material 
displacements u. With this obvious observation in mind, in (Bastrukov et al. 2009a, 2009b) 
focus was laid on non-investigated before regime of large lengthscale nodeless Alfven oscilla- 
tions, both global (in the entire spherical volume of star) and crustal (locked in the peripheral 
finite-depth spherical layer). The most conspicuous feature of this regime is that the radial 



^The decay time of equilibrium magnetic field of the neutron stars is much longer than the time intervals 
between X-ray bursts and periods of their quiescent pulsed emission (e.g. Bhattacharya & van den Heuvel 
1991, Chanmugham 1992, Goldreich and Reisenegger 1992), so that adopted approximation of an infinite 
electrical conductivity of the star matter is amply justified. 



3 



dependence of oscillating material displacements field u has no nodes. In a star undergoing 
global nodeless torsional oscillations, which are of particular interest for our present discus- 
sion, the fluctuating material displacements are described by the toroidal field of the form 
(Bastrukov et al 2007a, 2007b, 2009a) 

u(r, t) = [4>{r) X r]a{t), a{t) = aocosut, (5) 
5v{r, t) = u(r, t) = [6uj{r, t) x r], 6u){r, t) = 4>{r) a{t), (6) 

0(r) = Vx(r), V\(r) = 0, x{r) = j^r' Pe {cos 9) (7) 

where Pe is the Legender polynomial of degree i and the other symbols have their usual 
meaning. Fig.l illustrates the nodeless character of displacements in the star undergoing 
global non-radial differentially rotational, torsional, shear vibrations about polar axis in 
quadrupole and octupole overtones. 




Fig. 1. — Material displacements in the neutron star undergoing axisymmetric global torsional 
nodeless vibrations in quadrupole £ = 2 and octupole £ = 3 overtones. 

With the aid of the Rayleigh's energy method which is expounded below it was found 
that discrete frequencies of such vibrations is given by the following one-parametric spectral 
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formula (Bastrukov et al. 2009a) 



2M 

B _ 477 



1/2 



Va RB 



R V 3M 



^A = ^^, M=^pR' (9) 
V47rp 3 

the only parameter of which is the Alfven frequency, uja, of shear magneto-elastic oscil- 
lations of perfectly conducting stellar matter matter pervaded by magnetic field B in the 
star of radius R and mass M. It must be emphasized, however, that this theoretical spec- 
trum does not properly match the QPOs in the X-ray flux from flaring SGR 1806-20 and 
SGR 1900-1-14. One of reasons of this discrepancy may be inadequate assumption about 
homogeneous configuration of internal magnetic field and perhaps the most efficient way to 
clarify this conjecture is to investigate a model with geometrically different configuration of 
axisymmetric internal magnetic field. Before so doing it seems worth noting that the model 
of a star with homogeneous internal and dipolar external magnetic field has come into focus 
in astrophysics after seminal work of Chandrasekhar and Fermi (1953) in which the effect 
of mechanical flattening of the star at the poles of such magnetic field has been disclosed. 
Shortly after, similar conclusion has been drawn in outstanding paper of Ferraro (1954), 
but on the basis of star model with substantially nonhomogeneous internal and dipole-like 
axisymmetric poloidal magnetic field 

sm ti OU r sm u or 

U = Uin = ^r\?>r^ ~hR^)fiin^e, r < R, (11) 

U = Uex = —R^ sin^ 6, r> R (12) 
2r 

where B stands for the magnetic field intensity at the poles and V ■ B = as should be the 
cas^. The meridian cross section of Ferraro's model of the star is sketched in Fig. 2. 

Over the years, the different aspects, both mathematical and astrophysical, of Ferraro's 
model have been the subject of extensive investigations (e.g. Chandrasekhar & Prendergast 
1955, Robetrs 1955, Chandrasekhar 1956, Mestel 1956, Ledoux & Walraven 1958, Monaghan 
1965, Ledoux & Renson 1966, Sood & Trehan 1970, Goossens 1972, Goossens, Smeyers & 
Denis 1976). The effects of Ferraro's configuration of magnetic field (and magnetic fields 



^It may be noteworthy that magnetic energy stored in the star volume, W — (I/Stt) / dV, with this 
nonhomogeneous (nh) internal magnetic field Wnh — {69/252)B^R^ w 0.2AB^R'^ is somewhat larger than in 
the star with homogeneous (h) magnetic field Wh = {1/6)B^R^ w O.lTB^i?^. 



Fig. 2. — The meridional cross section of a neutron star with Ferraro's form of inhomogeneous 
poloidal internal and dipolar external magnetic field whose components are continuous on the star 
surface, contrary to a highly idealized star model with homogeneous internal and dipolar external 
magnetic field. 

of similar geometrical configuration) on the equilibrium shape, vibrational behavior and 
electromagnetic activity of pulsars and magnetars are considered in (Roberts 1981, loka 
2001, Braithwaite & Spruit 2006, Geppert & Rheinhardt 2006, Haskell et al 2008, Lee 2008; 
Broderick & Narayan 2008, see also references therein). 

In this work we focus on the non-studied before regime of global torsional Alfven nodeless 
vibrations of neutron star about axis of Ferraro's magnetic field f|T0l) - f|T2l) . In Section 2, 
the frequency spectrum of this toroidal mode is derived and compared with the frequency 
spectrum ([H]) of the neutron star model with homogeneous internal magnetic field. In Section 
3, the obtained spectral formula for the frequency is analyzed numerically in juxtaposition 
with data on QPOs during the flare of SGR 1806-20 and SGR 1900+14. Section 4 briefly 
accounts for the net outcome of this work. Technical details of analytic computations can 
be found in Appendix. 

2. Global Alfven torsional nodeless oscillations of neutron star in its own 

poloidal magnetic field of Ferraro's form 

In the model under consideration a neutron star is identified with a finite spherical mass 
of an elastic solid, regarded as an incompressible continuous medium of uniform density p 
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and an infinite electrical conductivity, whose vibrations under the action of Lorentz magnetic 
force are governed by equations of magneto-solid-mechanics ([2]) which can conveniently be 
represented in the following equivalent tensor form (e.g. Mestel 1999) 

pui = Vfc <5M,fc, SMiu = ]- \Bi 6Bf, + 5Bi - {B^ 5B^)5ik] , (13) 

47r 

55,(r, t) = {BkVk)ui - iukVk)B,, ViUi = (14) 

where 6Mik stands for the Maxwell's tensor of magnetic field stresses. The energy balance 
in the process of vibrations is controlled by equation 

J P^dV = - J SMikiiikdV = J [B, 5Bk + B^ 5B,] [Viiik + V^mJ dV, (15) 

u^ = UiUi, -itifc = -[VjMfc + VfcMj], Ukk = '^kUk = 0. (16) 

To compute the eigenfrequency of toroidal Alfven mode in question we take advantage of 
the Rayleigh's energy method which has been utilized in our previous above mentioned 
investigations. The key idea of this method consists in separable representation of fluctuating 
variables such as the vector field of material displacements Ui{r,t) and the tensor field of 
shear strains Uik{j,t) 

Ui{r,t) = ai{r)a{t), Uik{r,t) = ajfc(r) a{t), ajfc(r) = ^[Viafc(r) + Vfeai(r)]. (17) 

With this form of Ui, the magnetic flux density 6Bi{r,t) and the tensor field of fluctuating 
magnetic field stresses 6Mik{r,t) are represented in a similar manner 

6Bi{r, t) = 6,(r) a{t), 6,(r) = {BkVk)ai - {akVk)Bi, (18) 

6Mik{r,t) = [Tikir) - ir,/ifc]a(t), Tik{r) = -^[^.(r) 6fe(r) + ^^(r) bi{r)]. (19) 

The gist of this multiplicative decomposition of fluctuating variables is that on substituting 
(fT7|) - (fT9|) in (ITS!) this latter equation is reduced to equation for time-dependent amplitude 
a{t) having the well-familiar form 

Ma{t) + }Cma{t) = 0, M = jpaiaidV, (20) 
}Cm = J Tik fljfc dV = j [Bi bk + Bk bi] [V^ at + a^] dV. (21) 

Thus, from technical argument, the computation of frequency uj = [/C/TW]^/^ is reduced to 
calculation of integral parameters of inertia and stiffness JCm with the toroidal field of 
instantaneous, time-independent, displacements 

a, = A V X [rr'PeiO] : a, = 0, a, = 0, = A^il - (2)1/2^^1 (22) 
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and the magnetic field of Ferraro's form wliose splierical components inside tlie star are 

Br = l§-,r' (r^ - Ir'^ oos9, Be = (2r' - ^i?^) sin^^, B, = 0. (23) 




Fig. 3. — Frequency and period as functions of multipole degree i of global torsional Alfven 
vibrations of neutron stars with the Ferraro's shape of internal poloidal magnetic field. 

The torsional inertia as a function of multipole degree £ of nodeless differentially 
rotational vibrations in question is given by (Bastrukov et al, 2007, 2008) 



M 



+ 1) 



(2£+l)(2£ + 3)' 



(24) 



To avoid destructing attention from basic inferences of this work, we place all technical 
details of tedious but simple computations of integrals for /C^ in Appendix A. The final 
expression for this coefficient can be represented in the form 



kf 



(£2-1) (5£3 + 7f + 59£ + 84) 
2 (4^2 _ 1) + 3) {2i + 5) ' 



(25) 



And for the frequency spectrum of global nodeless torsional Alfven vibrations nodeless of 
the neutron star with Ferraro's form of nonhomogeneous internal magnetic field we obtain 



i^{oai) 



[i - 1) (5£3 + 7f + 59£ + 84) 



2 {2i - 1) {2£ + 5) 
VA B 



1/2 



va 



uja = B 



R 

3M' 



(26) 
(27) 
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It follows that the lowest overtone of this toroidal Alfven mode is of quadrupole degree, 
£ = 2. At £ = 1, the parameter of magneto-mechanical rigidity of neutron star matter 
cancels, /Cm(o«i) — 0' the mass parameter equals to the moment of inertia of rigid 
sphere, Ai{oa\) = J = {2/b)MR'^ . It follows from Hamiltonian of normal vibrations, 
H. = (l/2)7Wd^ + (l/2)/Ca^, that in this dipole case a star sets in rigid-body rotation, rather 
than vibrations, about axis of its dipole magnetic moment; this feature of the model under 
consideration is quite similar to that of the neutron star model with homogeneous internal 
magnetic field. In Fig 3., we plot the frequency z^(oa^) and the period P{oa\) = //^^(ga^) of 
the Alfven toroidal mode as functions of multipole degree computed in both homogeneous 
and nonhomogeneous neutron star models with indicated parameters. 

240 
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Fig. 4. — Theoretical fit (lines) of data (symbols) on low-frequency QPOs during the flare of 
SGRs 1806-20 and 1900-1-14 by the obtained spectral equation for the toroidal Alfven oscillations 
in Ferraro's poloidal field. 



3. QPOs in X-ray luminosity of flaring SGR 1806-20 and SGR 1900+14 from 

the viewpoint of considered model 

As was mentioned, the one-parametric spectral formula (IHl), computed in the neutron 
star model with homogeneous internal and dipolar external magnetic field, does not repro- 
duce general trends in data on QPOs frequencies whose numerical values for SGR 1806-20 
are given by z/jata =18, 26, 30, 92, 150, 625, 1840 and for the SGR 1900+14 these are 
'^data =28, 54, 84, 155 (Watts & Strohmayer 2007). It is tempting, therefore, to consider 
these data from the view point of investigated model by identifying the observed QPOs with 
overtones of spectral formula fl2^ . The result of i-po\e modal specification of the detected 
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QPOs frequencies as overtones of torsional Alfven seismic vibrations in question is presented 
in Fig.4 and Fig.5. Specifically, for SGR 1900+14 we obtain: u{oal) = 28 Hz; //(q^) = 53; 
Hz z/(oa*) = 84 Hz; 1/(0^13) = 155 Hz, and for the SGR 1806-20 we get z/(o4) = 26 Hz; 
z/(o4) = 30, z/(oa*io) = 92 Hz; z/(oa*6) = 155 Hz; //(ofl^g) = 625 Hz and z/(oa*8o) = 1840 Hz. 
While the lowest of detected oscillations, with u = 18 Hz, cannot be specified in terms of 
considered seismic vibrations, it is clearly seen that the obtained spectrum correctly reflects 
general trends in the detected QPO frequencies. This suggests, if the detected QPOs are 
really produced by Lorentz-force-driven global nodeless torsional seismic vibrations about 
the dipole magnetic moment of magnetars, their internal magnetic fields should be of sub- 
stantially nonhomogeneous configuration. 
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Fig. 5. — Theoretical fit (lines) of data (symbols) on high-frequency QPOs frequencies during the 
flare of SGRs 1806-20 by the obtained spectral equation for the frequency of toroidal Alfven mode. 

The result of £-pole modal specification of the detected QPOs frequencies as overtones of 
torsional Alfven seismic vibrations in question is presented in Fig.4 and Fig.5. Specifically, for 
SGR 1900+14 we obtain: z/(o4) = 28 Hz; z/(oa*) = 53; Hz z/(oa*) = 84 Hz; 1/(0^13) = 155 Hz, 
and for the SGR 1806-20 we get z/(o4) = 26 Hz; z/(o4) = 30, i^Mo) = 92 Hz; //(oalg) = 155 
Hz; 1/(0055) = 625 Hz and i^{oa\go) = 1840 Hz. While the lowest of detected oscillations, with 
u = 18 Hz, cannot be specified in terms of considered seismic vibrations, it is clearly seen that 
the obtained spectrum correctly reflects general trends in the detected QPO frequencies. This 
suggests, if the detected QPOs are really produced by Lorentz-force-driven global nodeless 
torsional seismic vibrations about the dipole magnetic moment of magnetars, their internal 
magnetic fields should be of substantially nonhomogeneous configuration. 
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4. Concluding remcirks 

Any attempt to predict the behavior of Alfven vibrational modes in pulsars and magne- 
tars, presuming of course that the star material possesses properties of perfectly conducting 
continuous medium pervaded by magnetic fields, is beset with uncertainties regarding ge- 
ometrical configuration of fossil internal magnetic field. It seems, therefore, that progress 
can be best made by studying these modes within the framework of comprehensive models. 
Among these are the models with homogeneous and nonhomogeneous axisymmetric poloidal 
magnetic fields considered long ago in works of Chandrasekhar and Fermi (1953) and by Fer- 
raro (1954), respectively, to show that such fields have the same effect as rigid rotation, that 
is, tend to produce a flattening of the star shape along the magnetic field axis. Following this 
line of argument and continuing investigations reported in (Bastrukov et al 2009a), we have 
computed here the frequency spectrum of axisymmetric torsional nodeless vibrations, in the 
neutron star model with Ferraro's form of nonhomogeneous poloidal magnetic field which is 
presented in Fig. 3 in juxtaposition in a neutron star model with homogeneous internal field. 

The practical usefulness of the obtained one-parametric spectral formula has been 
demonstrated by its application to £-pole identification of QPOs frequencies during the X-ray 
giant outbursts of SGR 1900-1-14 and SGR 1806-20. The result of our analysis, summarized 
in Fig.4 and Fig. 5, shows that the model adequately regains the overall trends in the detected 
QPOs frequencies and, thus, supports theoretical interpretation of these QPOs, advanced 
in works reporting this discovery (Israel et al 2005, Watts & Strohmayer 2006), as owing 
their origin to quake-induced torsional seismic vibrations of underlying magnetar. Together 
with this, in (Bastrukov et al 2009b) it has been shown that the same data on the QPO fre- 
quencies can be consistently interpreted from the view point of two-component, core-crust, 
model of quaking neutron star (Franco et al. 2000) with homogeneous internal magnetic 
field as being produced by axisymmetric differentially rotational Alfven nodeless oscillations 
of crustal solid-state plasma about axis of magnetic field frozen in the immobile core. With 
all these in mind, we conclude that it is the Lorentz force of magnetic field stresses plays 
decisive part in post-quake vibrational relaxation of above magnetars and that the toroidal 
fields of quake-induced material displacements are of substantially nodeless character. 

Authors are grateful to Dr. Dima Podgainy for helpful assistance. This work has been 
supported by NSC of Taiwan, grant numbers NSC-098-2811-M-007-009 and NSC-96-2628- 
M-007-012-MY3. 
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A. Torsional stiffness of global torisonal Alfven nodeless oscillations of a 
neutron star about axis of Ferraro's nonhomogeneous poloidal magnetic 

field 



In computing stiffness of torsional Alfven oscillations 

/Cm = y" Tik{r)aik{r)dV, aik{r) = ^[Vj afe(r) + aj(r)], 

T,fe(r) = ^[5^(r) bk{r) + (r)6,(r)], k{r) = (^^(r) Vfe)a,(r) - (afe(r) Vfe)5,(r) 
it is convenient to represent strain tensor 

in spherical polar coordinates with use of the angle variable C = cos^. In terms of this 
variable, the components of these tensor are 



dr 



aee = ttt H , 

r oQ r 



1 1 dttA, a,. 



c 



Ore 



r(l-C=^)i/2 Oct) r (1- (2)1/2 r ' 
1 
2 
1 
2 
1 



(1 — CY^"^ dttr ae dag 
r dC, r dr 
1 1 dar ^ da^ 



r (1 - C^)V2 d(t) r dr 
1 1 dae C (1 - ^y/^ da^ 



r (1 - (2)1/2 90 (1 _ (2)1/2 ^ r dC 

In the torsional mode of nodeless vibrations the field of instantaneous displacements has 
solely one non-zero (p — th component 

dm) 



a^ = ao = a^ = Ay{l-ey/'P^{C), P;(C) = 



In this case we have only two non-zero components of the strain tensor 



0-00 — are — 0, 



A 



(£-i)(i-C')V2p;, 



In spherical polar coordinates the components of vector field bi{r) = {BkVk)ai — {akVk)Bi 
are given by 
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Br 



dr r 

d__Bt 
dr r 

d ae 



d( r 



d_ 
84) 



de Be + B^ 



dC r 



ae + 



Bear-B^a^Ci'^-C)-^/'' 



-1/2 



d 



dr r 



dC r 



d(f) 



Be 



d ae 



_ _ (2)1/2 



-1/2 



dr r 



dC r 



d_ 
d^ 



Bj, — 



B^ar + B^ aeC(l - C 
r 

a^Br + a^BeCil-C)-'^^ 



Taking into account that Ferraro's field has only two non-zero components which can be 
conveniently represented in the form 



Br 



2/ 



C, Be 



(1 - CV^' 



f, B. = 0, f 



r^ r 
for the components of 6j we obtain 

br^O, be^O, b,^Br^-^{l-ey/' 

dr r 



da^ _ a^ Br _ a^BeCjl - C) 
d( r T 



where 



da^ 
dr 



Alr'~\l-eY''Pi 



da^ 



The integrand of Km reads 



'TikO'ik — 2 {Trfftartf, + Tecjittecj)) 

so that relevant to computation of /C^ components of tensor Tjfe = (l/47r)[Si bk + B^ 6i] are 
given by 



n4> 



1 

47r 
A 
An 
1 

At 



.BrBea^-ey 



da^ BrBe, >2xi/2^«</' 

^'^'^-^^^-^^ dc' 

{4{i - i)fV-'C{i - cY'Pe + 2//v-\(i - cY'i'^cp^ - + m]} , 

da^ BeBe , >2\i/2^^</' a^BeBr a^f, Be Be ({1 — Cy^'^ 

BeBr^ (i - C j -^7 

dr r dQ r r 



= ^ [-2/-^ {I - 1) //'c (1 - c') p; - r'-'r (1 - c') [2cp; - + 1) p^w . 
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The integral for stiffness can be conveniently represented in the form 

/C^ = 2 y [Tr^ ar4, + re^ a,^]dV = ^ {4(£ - ifRffh + 4(£ - l)Rff'[2Ii - i{£ + 1)12] 

+ Rf'f'i^h - M(£ + 1)72 + f(£ + Ifh]} 

The integrals Jj are computed with aid of standard recurrence relations between Legendre 
polynomials (e.g. Abramowitz & Stegan 1964) which yield 



J ai-e)p,p',dc 



2^(£+l)(2£^ + 2£-3) 
(4£2_i)(2£ + 3) ■ 

2i{i+l) 



(4£2-l)(2£ + 3)' 



/(I -a 



(4£2-l)(2£ + 3) 



For integrals with function / = (5/47?^) [7^2 ^3^,2 _ 5ii;2^j obtain 



rR r2 p2£+l 



Pff 



16 

^02^+1 

16" 



25 



30 9 

+ 



2£ + 1 2£ + 3 2£ + 5 



--R2, R2 — 



50 



90 36 
+ 



dr 



16 



2£ + 1 2£ + 3 2£ + 5 
100 



240 144 
+ 



2^+1 2£ + 3 2^ + 5 



The resultant expression for the stiffness reads 
where 



2 (4£2 - 1) (2£ + 3) (2£ + 5) 



(5^^ + 7^^ + 59£ + 84) . 
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